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Abstract. An improved generator for continuous unitary transformations is introduced 
to describe systems with unstable quasiparticles. Its general properties are derived and 
discussed. To illustrate this approach we investigate the asymmetric antiferromagnetic spin- 
1/2 Heisenberg ladder, which allows for spontaneous triplon decay. We present results for the 
low energy spectrum and the momentum resolved spectral density of this system. In particular, 
we show the resonance behavior of the decaying triplon explicitly. 
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1. Introduction 

Most low-energy properties of condensed matter systems can be understood in terms of 
quasiparticles, which are the elementary excitations of a system. Higher lying excitations 
are described in terms of scattering states or bound states of the elementary quasiparticles. 
One of the most successful concepts of quasiparticles was developed by Landau in the late 
1950's for interacting fermionic systems, the so-called Fermi liquid theory [1J. It is based 
on the idea that complex interacting excitations are adiabatically linked to the non-interacting 
ones. They share the same quantum numbers, i.e. momentum and spin. 

Only a few years later, Pitaevskii predicted that these quasiparticles can become unstable 
if certain decay channels exist [2J. The quasiparticles do not survive beyond a certain 
threshold in momentum space and their spectrum terminates at this threshold (see figured]). 
This prediction was confirmed later by neutron scattering measurements on superfluid 4 He 
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Figure 1. States with two excitations lying energetically below the single particle dispersion 
for momentum K < Kq- If in addition the Hamiltonian contains matrix elements which 
connect the one-particle space with the two-particle space, the (quasi)particles will become 
unstable for K < Kq (dashed line). The generator Fpc(l) leads to the dispersion consisting of 
the solid and the dashed-dotted lines (cf. section lXTT i. 

mm. 

In 2006 such a quasiparticle breakdown was measured for the first time in a quantum 
magnet [[H O. The spin excitations in the two-dimensional spin- 1/2 quantum magnet 
piperazinium hexachlorodicuprate (PHCC) show remarkable similarities with the excitations 
in superfluid 4 He. Stone et al. observed a threshold momentum beyond which the 
quasiparticle merges with the two-quasiparticle continuum and ceases to exist as well- 
defined excitation 0. This phenomenon was also observed in the quasi-one-dimensional 
antiferromagnet IPA-CUCI3 by Masuda et al. [6]. In theoretical considerations, magnons 
decaying into pairs of magnons are found in the long-range ordered Heisenberg model on 
the triangular lattice flTJ [H 13 as well as in the Heisenberg model on square lattices in strong 
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magnetic fields IfTOl . 

Generally, physical systems of unstable quasiparticles are much more common than 
systems of completely stable quasiparticles. For instance, in the case of Fermi liquids the 
quasiparticles always have a finite lifetime except at the Fermi energy. For the theoretical 
description and the understanding of many systems in condensed matter physics it is therefore 
an indispensable task to develop methods which are able to describe systems with unstable 
quasiparticles. 

A perturbative analysis of quasiparticle breakdown in quantum magnets was given in 
2006 by Kolezhuk and Sachdev IfTTTl and by Zhitomirsky lfT2l based on fully diagrammatic 
approaches. Both papers show that elementary excitations in gapped spin systems become 
unstable if they merge with the two-particle continuum. Spin systems in one and higher 
dimension are analyzed to explain the observations in IPA-CUCI3 and in PHCC. In one 
dimension a square-root dependence of the inverse quasiparticle lifetime is predicted [p~2|] . 
For the special case of an asymmetric rung-dimerized spin ladder Bibikov lfT3l confirmed 
these results by Bethe ansatz. An alternative approach to derive the lifetime of an excitation 
based on renormalization group methods was developed by Bach et al. lfT4l . 

Here we introduce an advancement of the method of continuous unitary transformations 
(CUTs) introduced in 1994 by Wegner [fT5l and independently by Glazek and Wilson lfT6l[T7ll . 
which allows us to describe systems with quasiparticle decay. 

The paper is divided into two main parts. In the first part in section El we give a 
short introduction to the method of CUTs and describe generally how one can deal with 
quasiparticle decay within this framework. In the second part consisting of section [3] 
and section HI we illustrate the general concept by explicit results for the asymmetric 
antiferromagnetic spin- 1/2 Heisenberg ladder. 

2. Method: Continuous unitary transformations 

In this section, we first outline the general concept of CUTs. Then we discuss similarities and 
differences between various schemes of CUTs depending on the choice of the generator. 

In principle, any Hamiltonian can be diagonalized by a suitable unitary transformation U. 
Famous examples are bosonization lfT8l or Bogoliubov transformations lfT9l whose fermionic 
version is used in the BCS theory of superconductivity EOl . For complex problems it 
is usually a very hard task to find such a suitable transformation. In 1994 Wegner [fT5l 
(and independently Glazek and Wilson US (13) presented a method to diagonalize a given 
Hamiltonian H in a continuous way. This method of CUTs is based on the idea to introduce a 
continuous auxiliary variable I and to define an /-dependent Hamiltonian H(l) := U\l)HU(l). 
Then the Hamiltonian transforms according to the flow equation 

d,H(l) = [F(l),H(l)], (2.1) 

with an anti-Hermitian generator F(l) := -U\l) (diU(l)). For / — » 00 the flow equation (12.11) 
maps the initial Hamiltonian H(0) := H to an effective Hamiltonian H e ^ := H(co) in a unitary 
way. Certainly, the final structure of the effective Hamiltonian // e g- depends on the form of 
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the chosen generator F(l). So the crucial point is to choose a generator F(l) which leads to a 
simplification of the initial Hamiltonian. Another important issue is whether the ensuing flow 
equation (12.11) is practically tractable. 

Wegner proposed to define the generator F(l) as the commutator between the diagonal 
part of the Hamiltonian H^(l) and the Hamiltonian H(l) itself. So the generator reads 
F(l) = [//^(Z), #(/)]. It was proven [fT51 l2TTl that this choice transforms the Hamiltonian in 
such a manner that [//^(oo), //(oo)j = 0, which implies that the final Hamiltonian H(oo) is 
block-diagonal with respect to the eigensubspaces of //^(oo). If H^(oo) is non-degenerate the 
final Hamiltonian H(oo) is actually diagonal. 

For band-diagonal Hamiltonian matrices, Mielke proposed another generator. His choice 
conserves the initial band structure during the flow 11221 . which is not the case for Wegner's 
generator. Mielke achieved the conservation of the band structure by introducing a sign 
function depending on the difference between the row index and the column index of the 
considered matrix element. 

Independently thereof, Knetter and Unrig [|23ll24l suggested a generator which allows us 
to create (quasi)particle number conserving effective many -body Hamiltonians. Their choice 
is also based on the idea to use a sign function. In contrast to Mielke's choice they used the 
difference of the particle number as the argument of the sign function. This generator can be 
regarded as a generalization of Mielke's generator for Hamiltonians formulated in second 
quantization. In the following, we denote this generator creating (quasi)particle number 
conserving effective Hamiltonians by Fp C (Z). An analogous generator was also used by 
Stein Il2~5ll26l for models where the use of the sign function was not necessary. 

In the following, we first summarize some properties of the generator Fp C (l) and specify 
its pros and cons. Particularly, we describe the problems arising in the description of systems 
with unstable (quasi)particles. Thereafter we present possible variations of the generator 
Fpc(/) including a generator which allows for the description of (quasi)particles with finite 
lifetime. 

2.1. The generator Fp C (l) 

Generally, a Hamiltonian in second quantization can be written as 



where H\(t) stands for the sum over all normal ordered terms which create i and annihilate 
j (quasi)particles, e.g. H^Q) is proportional to the identity and describes the vacuum 
energy during the flow. By the expression "term", we refer to both, the operators and 
the corresponding prefactor. The whole /-dependence of the Hamiltonian is carried by 
the prefactors. Note that for infinitely large systems the maximum number of involved 
quasiparticles N may be infinite, but this does not need to be the case. 



N 




(2.2) 



Adapted CUT to treat systems with quasiparticles of finite lifetime 



5 



According to the form of the Hamiltonian (12.21 ) the generator Fp C (/) is given by 

N 

F pc (l) = Y J Sgn(i-j)H i j (l). (2.3) 

i,j=0 

This means that terms in H(l) which contain more creation operators than annihilation 
operators are taken over to Fp C (l) with the same sign. Terms with more annihilation operators 
than creation operators are included in Fp C Q) with a negative sign. Terms leaving the number 
of particles unchanged do not occur in FpcQ)- 



For the generator Fp C (l) the flow equation (12.11) exhibits the following properties: 

a) If the spectrum of H is bounded from below, the flow equation converges |22ll24l| . This 
is the generic situation for physical systems. The mathematical derivation requires the 
Hilbert space of the system to be finite dimensional. 

b) The effective Hamiltonian H Q g is block-diagonal in the sense that it conserves the 
(quasi)particle number [|24l|. Therefore, the effective Hamiltonian commutes with the 
operator Q which counts the number of (quasi)particles 

[ffeff. Q] = 0- (2- 4 ) 
Thus it is of the form 



N 



tf e ff = 2^ (CX3) - (2 - 5 ) 

This property allows us to analyze subspaces with different (quasi)particle numbers 
separately. 

c) If the initial Hamiltonian H(0) has a block band-diagonal structure (i.e., Z/j(0) = for 
\i - j\ > N ), this block band-diagonal structure is conserved during the flow G21I241 . 

d) The generator Fp C (/) sorts the eigenvalues in ascending order of the particle number of 
the corresponding eigenvectors [|22ll27ll if the eigenvectors are linked by a matrix element 
of the Hamiltonian (see also section Appendix A). 

Items b) and c) are schematically illustrated in figure [21 Despite all the favorable 
properties of the generator Fpc^X it is not advantageous in every situation. Particularly, 
the last point is both a blessing and a curse. On the one hand, it ensures that the ground 
state is represented by the vacuum state of the effective mode|||. Additionally, it produces 
the appropriate (quasi)particle picture in systems where the elementary excitations have an 
infinite lifetime. But on the other hand, the described ordering of the eigenstates does not 
reflect the situation in many physical systems, e.g. systems with unstable (quasi)particles. 
This is schematically illustrated in figured] 

The generator Fpc(/) interprets the energetically lowest states above the ground state as 
the elementary excitations. In principle, it is possible to define the elementary excitations of 
the system in this way. But this definition can be misleading in the sense that states with very 
low or zero spectral weight are regarded as the elementary excitations of the system. Without 

$ For simplicity, the ground state |0> is assumed to be unique. 
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Figure 2. Schematical representation of the structure of the Hamiltonian H(l) during the flow 
for various generators. A colored block described by the pair i, j stands for the part H'.([) 
of the Hamiltonian. Only those blocks are colored where at least one term of H'Al) has a 
non-vanishing coefficient. In all cases, we assume an initial Hamiltonian which creates or 
annihilates at most two particles. For simplicity, we restricted our illustrations to terms which 
at most create or annihilate four particles for / > 0. Of course, terms which create or annihilate 
more than four terms may also occur. Panel (a) shows that the generator Fpc(/) conserves the 
block band-diagonality of the initial Hamiltonian during the flow and leads to a (quasi)particle 
number conserving effective Hamiltonian. Panels (b) and (c) show that both generators, Fgs(/) 
and Fg S do not conserve the block band-diagonality. The generator Fgs(/) only separates 
the H^(l) part, whereas the generator Fg S ip(0 also decouples the H\{1) part. 
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spectral weight we consider such states to be meaningless in terms of elementary excitations 
which serve as building blocks of all other excitations. Therefore, one usually defines the 
states with the largest spectral weight above the ground state as the elementary excitations of 
the system. Moreover, previous calculations Il28l |29ll strongly suggest that the rearrangement 
of the Hilbert space causes convergence problems in practice. In the perturbative approach of 
CUT E4ll30ll (p-CUT) these problems become perceivable in the extrapolations lf3Tll32~l . 

The second property, [# e ff, <2] = 0, of the effective Hamiltonian generated by Fp C (l) 
makes the describtion of unstable (quasi)particles difficult. By construction, the generator 
Fp C (Z) produces an effective Hamiltonian where the elementary excitations exhibit an infinite 
lifetime. The information of the decay is stored in the unitary transformation and therefore an 
additional transformation of observables is indispensable to describe the quasiparticle decay. 
This approach was first used by Kehrein and Mielke to describe dissipative quantum systems 

In the following subsection, we present a generator which does not eliminate the decay 
processes. Therefore, it is possible to study the quasiparticle decay more easily and more 
directly. The transformation of the observable is still necessary for quantitative results, but 
the essential aspect, i.e. the finite life time, is obvious without this transformation. 



2.2. Generator for the ground state 

To tackle the problems of (quasi)particle decay within the framework of CUTs mentioned in 
the previous section we introduce the adapted generator 

N 

F gs (D = Y J (H i (l)-H°(D) (2.6) 

relying on the form of the Hamiltonian (12.21) . We included only those terms in the generator 
Fgs(/) which either contain only creation operators or contain only annihilation operators. 
The terms which contain only creation operators are included as they appear in H(l). The 
terms which contain only annihilation operators are included with a negative sign relative to 
their sign in H(l). 

Again, the flow equation (12.11) converges if the spectrum is bounded from below. This 
follows directly from introducing a basis {\i)}, including the vacuum state |0), and examining 

0zffo,o(Q = -2£|fliv(9| 2 (2-7) 

m 

with Hij(l) := (i\H(l) \j). Note that Hij(l) describes an explicit matrix element in contrast 
to the previously appearing quantity H'Xl), which stands for a sum over terms in second 
quantization. According to (12.71) . H Q q(1) is a monotonically decreasing function of /. 
Therefore, if the spectrum is bounded from below, its derivative must vanish in the limit 
/ — > <x>. This also implies that 

lim H 0J (l) = lim H* (l) = 0, (2.8) 

i.e., all matrix elements connected to the vacuum state vanish in the limit I — * oo. In contrast 
to the generator i 7 pc(/) this generator destroys a block band-diagonal structure of the initial 
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Hamiltonian H(0). It solely separates the vacuum state from all other states. Hence the 
effective Hamiltonian is more difficult to analyze. This is the consequence of the more 
complex physics we have to describe. The evolution of the Hamiltonian H(l) during the 
flow using the generator Fg$(l) is compared to the one induced by Fp C (l) in figure [21 

While the choice (12.61) is very plausible, we have not presented a systematic derivation 
of Fg S (/) so far. To provide such an induction we adopt the derivation of a generator in 
the context of variational calculations [35 J. The idea of Dawson et al. was to minimize 
diE (J) = di (0| H(l) |0) under the constraint of a bounded F(l) so that the quantity E (l) 
decreases as fast as possible^. This leads to the calculation of 

<j{ (0| [F(0, H(l)] 10) + A ||F(/)||^ } = (2.9) 

with the Lagrange multiplier A > and denoting the Hilbert-Schmidt norm. With respect 
to a basis {\i)}, including |0), one obtains the expression 

4 Z ( F oADHi,o(l) - fliy(Ono(0) + * J] FIJI) F iV (/)| = (2.10) 

-FjAD 

with the matrix elements H t Jl) := (i\ H(l) \ j) and Fiji) := (i\ F(l) \ j). The variation implies 
= SojHjjo® - H Q ,i(l)ho ~ 2AF M (D (2.11) 

and hence 

F u (l) = ^ (HMSoj - S uo HoJD) ■ (2.12) 

In the following, we set A = 1 /2 and denote this generator by F m g S (/). It has the property that 
only matrix elements involving the vacuum state |0), i.e. Fiji) or F i (l), are different from 
zero. All other matrix elements vanish. 

The appealing property of i 7 m g S (/) is that there is a strong similarity to Fg S (t) in the 
sense that the terms of Fg$(l) containing only creation operators (or annihilation operators) 
represent the matrix elements Fiji) (or Fo,,-(Z)) of i 7 mgs(0 among other processes. But the 
effect on the total Hilbert space is very different. The matrix F m g S (/) is active if and only 
if there is a direct connection to the vacuum state |0), while, for instance, a term consisting 
only of creation operators also acts on states which already have a certain number of particles. 
Therefore, Fg$(l) can be seen as a generalization of i 7 mgs(0 for problems formulated in 
second quantization. But Fg S (Z) and i 7 m g S (Z) are not identical. 

The question arises if it is possible to adapt the above variational derivation of the 
generator i 7 mgs(0 to the generator Fgs(l) formulated in second quantization. This can be 
achieved by modifying the applied scalar product as we show next. 

We consider a system formulated in second quantization. Each operator acting on the 
Hilbert space can be represented by a sum over terms consisting of a product of creation 
and annihilation operators and a prefactor. We call the product of creation and annihilation 
operators a monomial. Thus a term consists of a monomial and a prefactor. 

§ To correspond with our approach in second quantization we use the vacuum state |0) as the starting vector for 
the minimization. In principle, one can use an arbitrary starting vector. 
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To obtain a unique representation of each monomial we first assume them to be normal 
ordered. Second, a certain ordering within all creation (annihilation) operators is implied. The 
creation and annihilation operators are denoted by e\ and e { , where 4 contains all quantum 
numbers describing the considered operator, for instance its position and spin. Note that 
such an expansion of a general operator is unique since all possible (ordered) monomials 
are linearly independent. They can be distinguished from one another by appropriate matrix 
elements. 

Next we define the scalar product of two monomials M\ and M 2 by 

( 1 for Mi = M 2 
(M U M 2 ):=\ (2.13) 
[ for Mi * M 2 

Since any operator on the total Hilbert space can be expanded in monomials, (12.131) in 
combination with the usual bilinearity of scalar products defines a valid scalar product. The 
scalar product (|2.13l) defines different monomials as pairwise orthogonal. So the set of all 
possible monomials are an orthonormal basis of the super Hilbert space of operators. 

The scalar product (12.131) implies the norm of an operator O as ||0|| 2 := (O, O). We again 
minimize (0| [F(t),H(l)] ]0), but with the constraint ||F(/)|| 2 = const. Thus we calculate the 
variation 

tf{ (0| [F(l), H(Q] |0) + A \\F(l)\\ 2 } = 0. (2.14) 
The operators H(l) and F(l) are expanded in second quantization 

H(l) = Yjh\(l)M\ (2.15a) 

and 

F{l) = Y J flH)M\ (2.156) 

(Ij) 

with the /-dependent prefactors [h]Q)) and {^'(O}- Here the bold indices i and j are sets of 
indices, e.g. i = fj'i, . . . , i N .}. Upper indices stand for creation operators and lower indices for 
annihilation operators. So Mj is short hand for the monomial 

M! = e) - --el e, ■■■e. . (2.16) 

The sums in (I2.15al) and (12.1561) run over all possible ordered sets i and j so that a unique 
expansion in monomials M! is achieved. 

Based on (|2.15al) and (12.1561) the right hand side of (|2.14l) to be varied has two additive 
contributions. The first one reads 

(0| [F(0, H(l)] |0) = (0| F(l)H(l) - H(l)F(l) |0> 

{i} 

where the empty set stands for the lack of non-trivial operators, in particular, a prefactor 
/; (/) belongs to a term that only contains annihilation operators. We exploit the fact that only 
creation operators yield non- vanishing results if applied to |0). Conversely, only annihilation 
operators yield non- vanishing bra states if placed right to (0|. 
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The second contribution reads 

^ ' (2.18) 



A(\\F(l)\\ 2 ) = A Yj\fl(l)\ 



Vi>J) 



Making the variation with respect to f!(l) vanish leads to 

Jj(f) = Yx - • (2.19) 

This generator solely contains monomials which are only composed of creation operators 
or only of annihilation operators. If we set A = 1/2 we obtain exactly the generator Fg$(l) 
we conjectured in (12.61) . Note that the above derivation holds for all kinds of operators in 
second quantization, including bosons, hard-core bosons, fermions and hard-core fermions. 
This terminates the derivation of Fgs(l) and its properties. 

In this paper, we only consider the case where the generator Fg S (Z) separates only the 
vacuum state from all other states. But we want to mention that it is also possible to generalize 
the generator Fg$(l) to the case where the vacuum state |0) is replaced by a statistical operator 
which defines a certain subspace, i.e. a reference ensemble. In this case the generator Fg$(l) 
induces an effective model on the reference subspace, which is separated from all other states. 
A well-known example is the derivation of the Heisenberg model or the t-J model from the 
Hubbard model. This generalization works very much in the same way as it was done for the 
generator F pc (/) before Il28ll29ll36ll37l. 

2.3. Other similar generators 

Besides the two choices of a generator considered so far (Fpc(/) in (12.31) and Fg$(l) in (12.61) ) 
there also exist other possibilities. For example, one can also include all terms to the generator 
FpcQ) which are connected to the one-particle subspace 

N N 

F gs,i P ( z ) = Z - + E (^ (/) - H ' (l) ) ■ (2 - 20) 

i>0 i>1 

Since this generator also separates the one-particle subspace from all subspaces with two and 
more particles, it is not an ideal choice to describe (quasi)particle decay. It suffers from the 
same caveats as Fpc(/). But this generator can be the optimal choice if the (quasi)particles 
have an infinite lifetime, while the higher particle subspaces are overlapping in energy (cf. 
figure [3]). In figure [2] the structure of the corresponding Hamiltonian H(l) is schematically 
illustrated during the flow. 

2.4. Common properties 

Although different generators produce different CUTs and therefore lead to different effective 
models it happens that they transform certain subspaces in exactly the same way. For example, 



it can be proven (see section | Appendix B.l[) that all generators considered in section 12.11 



section 12721 and section [2731 transform the vacuum state |0(/)) equally. This is a consequence 
of the fact that for all these generators the matrix elements from and to the ground state F i (l) 
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Figure 3. Overlap of the two- and three-particle continua. 



and Foj(l), respectively, are defined in the same way as long as the flow equation is treated 
exactly without any truncation. 

In figure |H we show numerical data verifying the equivalent transformation of the 
vacuum state |0(/)) by different generators. The /-dependence of the difference AE (l) : = 

1 

o.oi 

0.0001 

G- le-06 
< 

le-08 
le-10 
le-12 

5 10 15 20 5 10 15 20 

IJ IJ 

Figure 4. Evolution of AEo(l) := |£o(0 - £exact| during the flow for an antiferromagnetic 
spin-1/2 Heisenberg chain with 10 spins and exchange coupling J — J ± and x — 0, y — 1 
(cf. O.lal i- dTk/t ). All calculations started from the dimerized phase. The left panel shows the 
results for periodic boundary conditions; the right panel shows the results for open boundary 
conditions. 

\Eq(1) - £exact| between the vacuum expectation value E (l) := (0\H([)\0) and the exact 
ground state energy £exact ls plotted for the different generators. The system under study is 
an antiferromagnetic spin-1/2 Heisenberg chain with 10 spins and exchange coupling J. The 
starting point for all calculations is the ground state and the local triplons of the completely 
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dimerized phase (cf. section We considered periodic and open boundary conditions. 
Figure [4] shows clearly that all considered generators transform the vacuum state |0(/)) in 
the same way. The features beyond / « 12// ± stem from numerical inaccuracies occurring at 
AE « 10~ 10 /. These inaccuracies are shown here to illustrate where and how numerical errors 
make themselves felt. 

Similarly, one can prove that the generator Fp C (/) and the generator Fg S ^ transform all 
one-particle states identically (see section [Appendix ET2] ). 



3. Model: Asymmetric antiferromagnetic spin- 1/2 Heisenberg ladder 

The Hamiltonian for the asymmetric antiferromagnetic (J\\, J ± , /(jiag > ^) spin- 1/2 
Heisenberg ladder reads 

H = J ± (H ± + xH\\ + yH disLg ) (3.1a) 



with 



H ± = J] S ^ S 2,r (3.1*) 

r 

/ A (■Sl.r'S'l.r+l + ^>2,r^2,r+l) (3.1c) 
r 

/ A Si tr Sz,r+l, (3. Id) 



Hu 



^diag 

r 

where the first subscript 1,2 denotes the leg and r the rung (see figure [6]). The parameter x is 
given by x := J\\/J± and the parameter y by y := ^diag/Ti- 

This Hamiltonian contains some frequently discussed models. For example, for 
x = and y = 1 the Hamiltonian (|3.1al) - (l3.1JI) describes the exactly solvable isotropic 
antiferromagnetic spin- 1/2 Heisenberg chain JSC HS1 BD S3 H3 03- In the broad 
field of spin systems without magnetic long-range order the limit of the symmetric spin- 
1/2 Heisenberg ladder [|45l g6l SH HE [50l EH EH EH EH [551 HH EB (see figure © 
with y = is very popular as one-dimensional example of a valence-bond solid. Besides 
the theoretical interest in H in (I3.1al) - (I3.1JI) . there is a large number of compounds which 
can be described by spin ladders (see e.g. IfSTl [S8l [59l [60l f6Tl l62l l63l l64l 1651 l66ll or for an 
overview [|67l ). Special interest has been raised by the realization of coupled spin ladders in 
the stripe phases of cuprate superconducters [1681169117011 . Also the experimental evidence for 
superconductivity in Sro.ztCaB.eC^Ozn under pressure [71] contributed to the interest in the 
spin- 1/2 Heisenberg ladder and its extended versions. For the case x < y the Hamiltonian 
(13.1al) - (13.1<il) is usually denoted as dimerized and frustrated spin- 1/2 Heisenberg chain (see 
117211 and references therein). 

In the following, we study the two parameter sets x = 0.5, y = and x = 0.5, y = 0.1 
in the thermodynamic limit as generic examples. Since in both cases the relation x > y is 
fulfilled we call the system an asymmetric ladder instead of a dimerized and frustrated chain. 
A more comprehensive investigation of the dependences of the quasiparticle decay on the 
model parameter is left to ongoing research. The scope of this section and the next one is to 
illustrate the general considerations concerning the method by a concrete example. 
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Figure 5. Schematic representation of a spin ladder. Circles indicate spins with S = 1/2. 
Solid lines stand for couplings. The dashed line indicates the axis of reflection symmetry. 




r-2 r-1 r r+1 r+2 r+3 



Figure 6. Schematic representation of the asymmetric spin ladder. The additional diagonal 
interaction breaks the reflection symmetry and hence induces a hybridization between 

the one-quasiparticle states and the two-quasiparticle continuum. Thus it is an ideal model to 
study quasiparticle breakdown. 



The low energy spectrum for x = 0.5 and y = is well studied by several methods (see 
e.g. [47J and [1511 ') including methods based on CUTs [|73l . Therefore, it is a perfect starting 
point to discuss the more sophisticated case with x = 0.5 and y = 0. 1 . The additional diagonal 
interaction y makes the whole situation conceptionally more difficult because it breaks a 
symmetry. While for y = the model is symmetric under reflection (see figure [5]) an arbitrary 
small value y £ breaks this reflection symmetry (see figure [6]). 

The crucial point is that the reflection symmetry of the symmetric ladder is responsible 
for the infinite lifetime of the triplons, which are the 5 = 1 elementary magnetic excitations 
of an antiferromagnetic system without long-range order [74]. Breaking this symmetry 
creates processes which enable the triplons to decay into two-triplon states. Therefore, the 
asymmetric spin- 1/2 Heisenberg ladder is an ideal model to analyze quasiparticles with finite 
lifetime and to illustrate our previous theoretical considerations concerning the choice of an 
adaptive generator quantitatively. 

To define an appropriate starting point for the CUTs we use the bond operator 
representation llTBirToTl . Each rung r of the ladder considered separately has a four dimensional 
Hilbert space. A possible eigenbasis of the local operator 5i r 5 2 , r is given by the singlet state 

|j>,:=4=(IU>-liT», (3.2a) 
V2 

and the three triplet states 

tl r \s) r :=\t x ) r = ^-(m)-\ll)) r (3.2b) 
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(3.2c) 
(3. 2d) 



Without any interactions along the ladder (x = 0, y = 0) the ground state of the system is the 
product state of the rung singlets 

|0>:=f]|j> r . (3.3) 

r 

This reference state shall be the vacuum state of the system. Excitations on a rung r are 
created by the local operators t\ r , fj,,- and t[ r . These operators create a triplet on the rung r 
and satisfy the hard-core boson commutation relations 

lv» 4J = § r,s 8 afi - t\ r t ar - 8 afi J] ' (3 - 4) 

V y j 

where t a>r (a = x,y,z) annihilate such a triplet. We consider all the excited states, which 
can be continuously connected to the local triplets, to be the elementary magnetic excitations. 
They are called triplons 11741 1311 . 

In the bond operator representation the Hamiltonian (I3.1al) - (l3.1(il) is given in the notation 
ofd2j)by 

H = J ± (flg + H\ +Hl+Hl+H^ + H 2 l + flj) (3.5a) 
with 



r 

H\ = J] CUr + [\ X ~ \y) J] Z {Oa^l + ir+lQ 
r < ' r a 



\ z H ' r a±6 



a,r+l 



\ ^ r a 

Hi = (Hlf 

H * = ~\ y Z Z £ ^ (V + Wi)) 



r afi,y 



(3.5b) 
(3.5c) 

(3.5d) 

(3.5e) 

(3.5/) 
(3.5*) 

(3.5ft) 



This representation of the Hamiltonian of the asymmetric antiferromagnetic spin- 1/2 
Heisenberg ladder is used as the starting point for the CUTs. For y = the terms (3.5* ) and 
(13.5/? I) vanish whereby decay processes of one triplon into two are prevented. 
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Evaluating the commutator [F(l),H(l)] appearing on the right hand side of the flow 
equation (12.11) generates also terms which do not appear in the initial Hamiltonian H(0) = H. 
These terms must be added to the Hamiltonian with a coefficient equal to zero at / = 0. 
They must also be considered in the generator. Then the commutator [F(l),H(l)] generates 
even more terms which have to be taken into account. For a finite-dimensional Hilbert space 
this procedure comes to an end because the maximal number of terms is restricted. Such 
unrestricted calculations were performed to compute the results for the ground-state energy 
of the finite Heisenberg chain presented in section [2741 

For large systems such an unrestricted approach is not possible due to the proliferating 
number of terms. Especially in the thermodynamic limit, one has to deal with an infinite 
number of terms. Hence it is not possible to obtain a closed set of differential equations. Thus 
in practice one has to decide which terms are important to describe the underlying physics 
properly and which terms can be neglected. 

One established truncation scheme is to use a perturbative approach 11241 1301 (p-CUT) 
which is based on the generator Fpc(/). But since we intend to describe the decay of 
quasiparticles so that variations of the generator Fp C (l) (see section |2]) have to be used, we 
choose the self-similar approach (s-CUT). But there is no fundamental reason why the adapted 
generator cannot be implemented perturbatively as well. 

The s-CUT was used in many previous applications of the CUTs (for an overview see 
[1771 and references therein) among them the original work by Wegner [fT5l . The whole 
transformation takes place in the coefficients of the terms in the Hamiltonian which motivates 
the naming "self- similar". This approach can straightforwardly be implemented for various 
generators. 

In the present paper, we apply the following truncation scheme, which is based on the 
finite correlation length of one-dimensional systems with a finite gap. The truncation scheme 
is based on the locality of the term which is justified for systems with finite correlation length. 
It is described in detail in the following section [3TT1 

All in all we proceed in the actual calculations as follows: 

• Define a truncation scheme which restricts the maximal number of terms. 

• Set up the flow equation by calculating the commutator [F(l), H(l)]. Only those terms 
are considered which fit the truncation scheme. 

• Solve the flow equation numerically. 

For calculations in the thermodynamic limit, it is necessary to make use of the 
translational invariance of the Hamiltonian (I3.5al) - (|3.5/zl) . The translation symmetry ensures 
that terms which describe identical processes except for a shift along the ladder have the same 
coefficient. Consequently, it is sufficient to track only one representative of this symmetry 
group. This procedure is also possible for all others symmetries of the Hamiltonian (I3.5al) - 
(13.5/? I) , e.g., the spin symmetry and the rotational symmetry by n of the ladder. Using 
representatives of the underlying symmetries reduces the number of coefficients appearing 
in the flow equation (|2.1I) significantly for a given truncation scheme. Thus more extended 
truncation schemes considering more processes become feasible. 
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3.1. Truncation scheme 

A truncation scheme is necessary to limit the number of terms so that a closed set of 
differential equations is achieved. In the present paper, we are dealing with terms in real-space 
and the truncation scheme is based on the locality of the terms. We first define a measure for 
the locality of a term, which we call the extension d. The extension d of a term is defined 
by the distance between the rightmost to the leftmost rung on which the monomial acts in a 
nontrivial way. For example, the term t\ r t^ r+1 t r+4 has an extension d = 4. Second, we define 
the truncation scheme by choosing a maximal extension rfmax discarding all terms with a 
larger extension (d > rfmax)- 

It turns out that it is appropriate to define not only one maximal extension for all terms 
but to keep terms with a different number of annihilation or creation operators up to different 
maximal extensions ||29l . Accordingly, terms with n annihilation or creation operators in 
total are required to have an extension d n or less to be kept in the flow equation. As a 
second truncation criterion we admit only terms which create or annihilate not more than 
N (quasi)particles. Thus the total truncation scheme is defined by the value of N and the 
set of extensions d = (J2, ■ • ■ , ^2a0- Note that due to the conservation of spin no single 
triplon operators can occur in the Hamiltonian. So no d\ needs to be denoted. In addition, 
the translational invariance of the Hamiltonian makes d\ superfluous. In this case only six 
different monomials exist which act on one rung only, in particular, 4,r and t a>r with a = x, y, z. 
For the symmetric ladder (y = 0) no terms occur which consist of an odd number of operators. 
Therefore, we do not need to define maximum extensions J 3 , J 5 , . . . . In the notation of the set 
of extensions d we replace such superfluous extensions by a dot, e.g. d = (8, ., 6, ., 4). 

It is worthwhile to emphasize that this truncation scheme does not turn our approach to 
a calculation on a finite cluster. It is a self-similar calculation strictly in the thermodynamic 
limit. We only truncate the range of the interactions in real space, but not the Hilbert space. 

4. Results 

Here we present results for the symmetric ladder with x = 0.5 and y = and for the 
asymmetric ladder with x = 0.5 and y = 0.1. These two parameter sets are chosen to 
illustrate the differences between systems with stable quasiparticles and systems with unstable 
quasiparticles which exhibit a finite lifetime. In particular, we confirm our previous statements 
concerning the properties of the different generators (cf. section [2]). 

Firstly we show that a rearrangement of the states of the Hilbert space, i.e., a continuous 
re-labelling (for simple examples see EU), reduces the speed of convergence (see section HTTT) . 
Therefore, generators which avoid such a rearrangement induce a considerably faster 
convergence. Secondly, we discuss the low energy spectrum for the symmetric and for the 
asymmetric ladder (see section 14.21) . If decay is possible the generator Fg S ip(Z) and the 
generator Fp C (/) indeed tend to interpret the energetically lowest states above the ground 
state as the elementary excitations (as stated before in section [2]). This can be avoided by 
using the generator Fg S (Z). Unfortunately, for this generator a simple calculation in the one- 
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particle subspace is not sufficient to arrive at reliable results for the true one-triplon dispersion. 
This is a consequence of the fact that in the effective Hamiltonian induced by Fg$(l) the one- 
particle subspace still couples to higher particle subspaces (see figure ©. To obtain reliable 
results for the single triplon dispersion we include states which consist of up to three particles 
in our calculations (see section |4~3T) . Especially, we show results for the zero temperature 
spectral density in which the (quasi)particle decay is manifest as a Lorentzian resonance of 
finite width. 



4.1. Convergence 

To quantify the speed of convergence of the flow equation for different generators we 
introduce the residual off-diagonality (ROD) Il28l |29ll . The ROD is defined as the square 
root of the sum of the moduli squared of all coefficients that contribute to the considered 
generator. Using the notation of (12.151? I) the ROD is given by 



ROD(/) 



(4.1) 



where the range of the sum depends on the choice of the generator. Note that only one 
representative of the translational symmetry group is included. Otherwise the ROD would 
grow proportional to the system size. In addition, ROD^ denotes the square root of the sum 
of the moduli squared of all coefficients belonging to terms with i creation and j annihilation 
operators or to their Hermitian conjugate terms. 

Figure shows the evolution of the ROD during the flow for different generators and 
different truncation schemes for x = 0.5 and y = 0. For all generators the RODs decrease 
strictly monotonically. The ROD of the generator Fg$(l) decreases faster than the ROD of the 
generator Fg S ip(/). This is a consequence of the fact that the generator Fg S ip(Z) contains 
more coefficients than the generator Fgs(/). The convergence of these additional coefficients 
is slower because they connect states which differ less in their eigenenergies (cf. (IA.6I) ). e.g. 
the energy gap between one- and three-triplon states is smaller than the energy gap between 
the vacuum state and the two-triplon states. This also explains why the generator Fg S ip(/) 
converges faster than the generator Fp C (/). 

The convergence behavior clearly changes if one includes the diagonal interaction, even 
if y is small (y = 0.1). In figure |7b the ROD during the flow for different generators and 
different truncation schemes for x = 0.5 and y = 0.1 is depicted. Only the ROD of the 
generator Fg S (/) decreases strictly monotonically while the RODs of the generator Fg S 
and the generator Fpc(/) increase temporarily during the flow. These increases indicate a 
rearrangement in the Hilbert space, cf. ETTl for simple examples. If all eigenstates were 
ordered in such a way that states with more triplons had higher eigenenergies the ROD would 
decrease exponentially (cf. (IA.5I) ). These rearrangements affect the results for the one-triplon 
dispersion as we illustrate in the section I4~2l 

Figure|7]: shows the ROD of the generator Fpc(/) split in the partial RODs ROD', defined 
above for x = 0.5 and y = 0.1. Clearly, the contributions ROD^ of the ROD changing the 
number of triplons only by one (\i - j\ = 1) provide the main contributions to the total ROD, 



Adapted CUT to treat systems with quasiparticles of finite lifetime 



18 




ROD, F ,N=4, d=( 10,8,8,5,5,3,3) 



ROD, F , N=5. d=(12,10,10,6,6,5,5,4,4) 
ROD, F , N=5, d=(18,14,14,6,6,5,5,4,4) 



100 
/■./ 



150 



200 



200 



Figure 7. Convergence of the flow equations. Panel (a) shows the evolution of the ROD 
during the flow for different generators and different truncation schemes for the symmetric 
ladder (x = 0.5 and y = 0). In all cases the ROD decreases strictly monotonically. Panel 
(b) shows the evolution of the ROD for different generators and different truncation schemes 
for the asymmetric ladder (x = 0.5 and y = 0.1). The RODs of the generator Fg S ip(0 
and the generator Fpc(l) increase temporarily during the flow. This indicates a significant 
rearrangement of the states in the Hilbert space. Panel (c) shows the ROD of the generator 
Fpc(l) split in the parts ROD^ for the asymmetric ladder (x = 0.5 and y = 0.1). The main 
contributions to the total ROD is due to ROD' with |/ - j\ — 1. 



although the corresponding initial couplings are proportional to y, which is small (y = 0.1). 
From this we infer that the convergence of the flow equation is mainly influenced by terms 
which induce a rearrangement of the Hilbert space if they are to be eliminated by the CUT. 
It is less important whether the corresponding coupling parameter is large or not. This is an 
important property of the CUTs which distinguishes them from conventional diagrammatic 
perturbation theories. 
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In summary, we state that a rearrangement of the states of the Hilbert space reduces the 
speed of convergence. Omitting the corresponding terms from the generator stabilizes the flow 
in the sense that the convergence is monotonic and robust. Hence, especially the generator 
Fgs(/) yields a very fast converging and robust flow. We point out that a fast convergence 
is advantageous because it minimizes the interval in / during which significant terms are 
truncated. Hence as a rule of thumb, the faster the convergence, the smaller are the truncation 
errors. 



4.2. Low energy spectrum 

Here we discuss the low energy spectrum of the effective Hamiltonian // e ff. We always 
stopped the CUT at / = 200/7 ± . At this value the remaining effect on the one-particle subspace 
is small in all cases (cf. figure [7]) so that a further integration of the flow equation would not 
change the results for the one-triplon dispersion as shown in figure [8] 

To calculate the one-triplon dispersion u>\(K) of the effective model we define the 
Fourier-transformed one-particle states 

\z> a )-=-j=Yj Q{Kr \ r > a) (4 - 2) 

with \r, a) := 4,r 10). The action of H\(l) with respect to the translational symmetry is given by 



(ID. ID in section [Appendix D~[ Due to the SU(2) symmetry of the Hamiltonian H(l) the hopping 
coefficients c a n ' a r obey the relation c a n ' a r = 6 a ^ a Cn- r . This leads to the threefold degenerate one- 
triplon dispersion 

oi X (K) := (K, a\ H{(oo) \K, a) = J] ^cn-r ■ (4.3) 

r 

The generator Fpc(/) and the generator Fg S ip(/) separate the one-particle subspace 
from the other subspaces. Consequently, for these two generators the one-triplon dispersion 
co\{K) yields eigenvalues of the effective Hamiltonian H Q ff in the one-particle subspace. In 
contrast, the generator Fg S (Z) does not separate the one-particle space. Therefore, the effective 
Hamiltonian H q q still contains terms which connect the one-particle subspace with higher 
particle states. In this case the quantity a>\(K) only gives an approximation of the eigenvalues 
of the effective Hamiltonian H e ^ (cf. figure IDTV ). 

In figure [8k the one-triplon dispersion a>i(K) is displayed for x = 0.5 and y = 0. 
Results for all three generators Fpc(/), Fg S ip(/) and Fg$(l) and various truncation schemes 
are shown. The two generators Fp C (/) and Fg S ip(Z) separating the one-particle space yield 
almost the same results and barely depend on the chosen truncation scheme. Together with the 
good convergence (cf. figure[7k) this implies that the results are very reliable. By construction, 
for the generator Fgs(/) the quantity u)i(K) as defined above yields only an approximation of 
the true one-triplon dispersion. The resulting o)i(K) is an upper bound to the results obtained 
from the other two generators if the truncation errors are negligible. This fact is based on 
the variational principle that a minimum in a restricted subspace is an upper bound to the 
minimum in an unrestricted subspace. To improve the results in this case one has to consider 
higher particle subspaces as well (cf. section |4~3~T) . 
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Figure 8. Low energy spectrum of the symmetric and asymmetric spin-1 /2 Heisenberg ladder. 
Panel (a) shows results for the one-triplon dispersion u>\(K) for the symmetric ladder with 
x = 0.5 and y = 0. Results for different generators and different truncation schemes are 
depicted. Additionally, the lower part of the two-particle continuum is shown (grey area). The 
dashed white line represents an approximation of the lower edge of the two-particle continuum 
obtained by the approximate one-triplon dispersion to\{K) in the case of the generator Fg$(l). 
Panel (b) shows the corresponding quantities for the asymmetric ladder with x = 0.5 and 
y = 0.1. 



Figure [8] also displays the lower part of the two-particle continuum 

co 2 (K, Q) = co, (1/2 + Q) + co l (K/2 - Q) , {A A) 

where Q e [-n, n] denotes the relative momentum. The lower band edge is given by 
the minimum of (x>i(K, Q) over Q; the maximum yields the upper band edge, respectively. 
Here we used the one-triplon dispersion ^i(^) we obtained by the generator Fp C (l). The 
additional dashed white line represents an approximation of the lower edge of the two- 
particle continuum obtained by the approximate one-triplon dispersion a>i(K) in the case of 
the generator Fgs(l). We emphasize again that due to the reflection symmetry for y = (cf. 
figure [5]) no interaction exists between the one-triplon states and the two-particle continuum. 
As a result the quasiparticles are well-defined and infinitely long-lived for the whole Brillouin 
zone, although the two-particle continuum starts below the one-triplon dispersion for certain 
momenta K. In addition, this symmetry prevents any rearrangement between the one- and 
two-particle subspaces during the flow (see section Appendix A). This situation changes 
abruptly if a diagonal interaction is switched on, even if y is very small. 

In figure [8h the one-triplon dispersion coi(K) is displayed for x = 0.5 and a small 
additional diagonal interaction y = 0.1. Again results for all three generators FpcO), Fg S 1 p(/) 
and Fg$(l) and various truncation schemes are shown as well as the lower part of the two- 
particle continuum a>2(K, Q) determined from the one-triplon dispersion obtained by the 
generator Fp C (l). Likewise, the approximate results for the lower edge of the two-particle 
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continuum obtained by the generator Fg S (/) are shown. For comparison only, we also included 
the former results of the one-triplon dispersion obtained by the generator Fpc(/) for x = 0.5 
and y = 0. 

The use of the two generators Fp C (l) an d Fg S ±^(1) implies significantly lower energies 
for the one-triplon dispersion, see figure [8b, where io\(K) overlaps with the two-triplon 
continuum. The results strongly depend on the truncation scheme in this region. This can 
be explained as follows. Since for y ± the one-particle and the two-particle space are 
interacting with each other, the two generators Fp C (/) and Fg S try to sort the eigenvalues 
in such a way that the eigenvalues of the one-triplon dispersion lie below the two-particle 
continuum, see figure [Hand section [Appendix A} Therefore, the one-triplon dispersion of the 



effective model lies at the lower edge of the two-particle continuum in the region where 
the one-triplon dispersion merges with the two-particle continuum. This is not completely 
achieved in practice because of the indispensable usage of a truncation scheme. 

We truncate the range of the decay processes in real space. This means that the distance 
between the generated two triplons is limited although the true scattering state comprises 
contributions up to infinite distances. As a result, the rearrangement of the eigenvalues is 
only incomplete. Figure [8b illustrates that increasing the range of the decay processes (e.g. 
increasing <i 3 ) implies that oj\(K) approaches the lower band edge of oj\(K, Q) from above 
more and more. 

As stated before, the rearrangements of the states are unfavorable for two reasons, (i) 
They imply a slow convergence which may cause growing truncation errors, (ii) One usually 
defines the state with the largest spectral weight above the ground state as the elementary 
excitation of the system and not a state with almost no spectral weight, even if it is lower in 
energy. 

To avoid the rearrangement of the eigenstates, which leads to a potentially misleading 
quasiparticle picture, we employ the operator Fg$(l) (cf. figure [8b). As before the 
generator Fg$(l) only yields an approximation for the one-triplon eigenvalues of the effective 
Hamiltonian H Q ^. This is the case even in the region of the Brillouin zone where the 
quasiparticles are well-defined. Due to our treatment of the problem in real space we cannot 
distinguish processes in different regions in momentum space easily. To improve the results 
for the one-triplon dispersion one must consider the interaction with states which consists of 
more than one particle as well. This is discussed in the next section 1431 

Here, we first want to show the results for the two- and three-particle continua resulting 
from the approximate one-triplon dispersion a>i(K) in the case of the generator Fg$(l) for 
x = 0.5 and y = 0.1. The three-particle continuum can be determined in analogy to the two- 
particle continuum. One only has to replace one one-triplon dispersion on the right hand side 



of (14.41) by energies of the two-particle continuum. 

The boundaries of these continua are shown in figure 9] by solid lines. Additionally, this 
figure shows the boundaries of the three-particle continual emerging from the combination 
of the approximate one-triplon dispersion a)\(K) and a certain approximate two-particle (anti- 

|| The name three-particle continuum refers to the fact that the corresponding states consist of three triplons in 
the basis of the effective Hamiltonian. 
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Figure 9. Two- and three-triplon continua of the asymmetric spin-1 /2 Heisenberg ladder with 
x = 0.5 and y = 0.1. The solid green lines represent the lower and upper boundaries of the 
two-particle continuum. The other green lines represent two two-particle bound states and 
one two-particle anti-bound state. The blue lines illustrate the boundaries of the three-particle 
continuum, where the type of the lines correspond to the two-particle state(s) used to determine 
the three-particle continuum. 



)bound state. Here we only use the region where the corresponding (anti-)bound state is 
well-defined, i.e., does not merge with the continuum. 

The two-particle bound states and the two-particle anti -bound state shown in figure [9] are 
calculated by diagonalizing the effective Hamiltonian H e ^ in the subspace spanned by the 
single-triplon states 

\K,a):=-}=y \^ r \r,a) (4.5a) 

and the two-triplon states 

\K, a) \d,J3) := V e 1 ^) \r, a) \r + d,j3) (4.5b) 



with < d < 120 for each given value of K (for details see section [Appendix D| 

Since the subspace spanned by the states (|4.5d) and (|4.5l?l) is not separated from higher 
triplon states (cf. figure ID lb we obtain - as for the one-triplon dispersion a>i(K) - only 

1 Here we considered d < 120 only to be consistent with the later calculations which also include the three- 
particle space. 
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an approximation for the bound states. Consequently, the depicted continua only represent 
approximations as well. The restriction of the relative distance d in the two-triplon subspace 
(14.51? I) is less important. Increasing d does not change the results perceivably. 

As in the symmetric case x = 0.5 and y = EH EH [801 S3 ED EH, two two-particle 
bound states - one with total spin 5=0 and one with total spin 5 = 1- and an anti-bonding 
state with total spin 5=2 exist. The boundaries of the continua shown in figure [9] help to 
understand the shape of the spectral densities calculated in the following section. 

Note that the whole complex structure of the low energy spectrum shown in figure [9] 
follows from the one-triplon dispersion, the triplon-triplon interaction, and from the 
diagonalization of the effective Hamiltonian within the subspace (14.5al) - (14.5&[) . 



4.3. Spectral density 

In this subsection we improve the results presented in the former section for the one-triplon 
dispersion which we obtained by the generator Fg S (Z) for x = 0.5 and y = 0. 1 . This is achieved 
by including interactions with three-triplon states. To describe triplon decay we calculate the 
zero temperature spectral density. 

We start by analyzing the frequency and momentum resolved retarded zero temperature 
Green function 

G (K, co) = Urn (K, z\ [co - (# eff - E ) + i6]~ l \K, z) . (4.6) 

The spectral density 5 (K, co) follows by taking the negative imaginary part of G (K, co) divided 
by 7r 

5 (K, co) = --3 [G (K, co)] . (4.7) 
n 

The Green function is evaluated by tridiagonalization (Lanczos algorithm) which leads to the 
continued fraction representation lf8Tl 1821 18311841 l85ll 



G(K, co) = 1 7 ■ (4.8) 

" — Sfw 

u-ai(K)-— 

The coefficients a t (K) and bi(K) are calculated by repeated application of H q q - E on the 
initial state \K,z) with wave vector K, spin 5 = 1, and 5, component m = (for details see 
section Appendix C). Note that the continued fraction in the denominator on the right hand 
side of (14.81) (proportional to bi(K) 2 ) can be taken as a standard self-energy whose imaginary 
part determines the decay rate. In this respect our approach is not so different from the one in 

ma. 

In all practical calculations we have to restrict ourselves to a certain subspace. For this 
calculations we considered the subspace spanned by 

\K,a):=-^-Ye iKr \r,a), (4.9a) 
\K, a) \d,/3) := J= £ e 1 ^) \r, a) \r + d,J3) (4.9b) 
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and 

1 V iK(r+ 2dl+d2 ) 

\K,a)\d u /3)\d 2 ,y) 2_j e [ ~^~ , \r,a)\r + dufi^\r + d 1 +d 2 ,y) (4.9c) 

with d, d\, di > and d, d\+d 2 < 120. Note that (I4.9al) is the one-triplon state for fixed K and 
a, (14.961 ) the two-triplon scattering states, and (|4.9cl ) the three-triplon scattering states. Thus 
we only need the restricted effectiveHamiltonian 

= J ± (H{(oo) + Hj(oo) + # 3 3 (oo) + H 2 (oo) + Hl(pa) + H\(oo) + H$(po) 

+ // 2 3 (oo) + // 3 2 (oo)) <4 ' 10> 

The action of this restricted effective Hamiltonian on the subspace (I4.9al) - (|4.9cl) is 



presented in section Appendix D More details about the calculation of the spectral density 



are given in section Appendix C 



Figure [10a shows the spectral density S(K,co) for x = 0.5 and y = 0.1. To keep track 
of this spectral density we display the results presented in figure [9] also in the co, K-p\ane of 
figure [10a. For K = and K = n/\0 decaying triplons are observed. Their density lies in 
the vicinity of the approximate one-triplon dispersion. The region framed in red is shown in 
detail in figure [10b. In this region we can fit our data by a Lorentzian 

A (F/2) 2 

Ufo) = 7 (4. 1 la) 

n(oj-oj ) 2 + (T/2) 2 

with 

A» 1.0848/ Jx (4.116) 

£1)0 » 1.4938 J ± (4.11c) 
and the inverse lifetime 

r«0.04027 ± . (4. lid) 

For clarity, figure [10b shows the spectral density at K = n/lO on logarithmic scale. 
Besides the strong one-triplon peak at o) « 1.49387 ± the spectral density increases distinctly 
at the beginning of the three-particle continuum involving the 5=0 bound state. Another 
rise of the spectral density occurs at the upper end of the three-particle continuum involving 
the 5 = 1 bound state. At the beginning of the three-particle continuum involving the 5=2 
anti -bound state the spectral density drops notably. This illustrates that the existence of bound 
and anti-bound states influences the form of the spectral density significantly. Note that an 
additional channel, for instance contributions from scattering states of a single triplon and a 
triplon-triplon bound state, does not always imply an increase of the spectral density. It may 
also lead to a significant decrease. We attribute this phenomenon to destructive interference. 
That means the additional channel interferes destructively with the already existing channel 
so that a net decrease is engendered. 

Finally, we want to discuss the shift of the one-triplon dispersion caused by the 
hybridization of two- and three-particle states (cf. (|4.9al) - (14.9cl) ). In the sequel we call this 
shifted one-triplon dispersion as the renormalized one-triplon dispersion. The results shown 
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Figure 10. Spectral properties of the asymmetric spin- 1/2 ladder. Panels (a), (b) and (c) show 
the spectral density of the asymmetric spin-1/2 ladder for x = 0.5 and y = 0.1. In Panel (b) 
a Lorentzian is fitted to the spectral density. The described quasiparticle exhibits a inverse 
lifetime of T « 0.0402/ ± . Panel (d) compares the results of figure [SJ) with the red crosses 
depicting the renormalized one-triplon dispersion obtained by using the generator Fgs(0 an d 
(tri)diagonalization in the subspace (14.9ab -( l4~9c] i. This renormalized one-triplon dispersion is 
also depicted in Panel (a). 

in figure [TOli are partly those of figure [8b- In addition, the red crosses depict results obtained 
from a tridiagonalization after the CUT induced by Fgs(l). 

In the region of the Brillouin zone where the quasiparticles are well-defined we obtain 
the renormalized one-triplon dispersion by fixing the total momentum K and calculating the 
lowest eigenvalue of the effective Hamiltonian in the subspace (14.9d) - (14.9cl) . The results 
are represented in figure [TOa and figure ITOti by red crosses. To good accuracy, we retrieve 
the results obtained before by Fpc(/) and Fg S ip(Z) in the region without decay. Therefore, 
it is sufficient to consider the subspace (|4.9al )- (|4.9cl ) if one wants to describe the one-triplon 
dispersion of the asymmetric spin-1/2 Heisenberg ladder with x = 0.5 and y = 0.1 by using 
the generator Fg S (/). 

Note that the calculation in the subspace (|4.9al )- (|4.9cl ) does not lead to the correct band 
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edges of the triplon continuum because the shift of the one-triplon dispersion makes itself felt 
only if we included four-particle states as well. Hence this kind of calculation is not fully self- 
consistent. There are possibilities to achieve consistency between the one-triplon dispersion 
and the band edges of the continua. But this issue is beyond the scope of the present article. 

In the region where the one-triplon dispersion hybridizes with the two-triplon continuum 
the red crosses shown in figure [10a and figure ITOU indicate the energy with the maximum 
spectral intensity S(K,to). These energies represent what is usually seen as the energy of a 
quasiparticle with finite life-time. The energies determined in this way lie between what is 
obtained from Fg$(l) in the one-triplon sector (blue line in figure [TObD and what is obtained 
from Fpc(/) or from ^g Sj ip(0- 

We emphasize, that the advantage of the generator Fg$(l) compared to the generator 
Fg S ip(/) or Fpc(/) is that also the quasiparticle decay is described in the region of the 
Brillouin zone where the one-triplon dispersion merges with the two-triplon continuum. The 
generator Fg S (l) avoids rearrangement processes during the flow which lead to a potentially 
misleading quasiparticle picture. Thereby the CUT becomes more robust. Hence the proposed 
adapted generator achieves the goal from the outset to describe decaying quasiparticles 
properly. 

5. Summary 

In the present paper, we introduced an approach based on continuous unitary transformations 
to describe systems with unstable quasiparticles. The main idea is to use a generator 
formulated in second quantization which leads to an effective model where only the ground 
state is isolated from the rest of the Hilbert space. 

In the first part of this paper we described the properties of this adapted generator 
and discussed similarities and differences with other generators. Additionally, we derived 
the adapted generator in the context of variational calculations. All considerations were 
completely general and did neither depend on the model under study nor on the actual 
realization of the continuous unitary transformation. Thus we expect that generally an 
analogous modification of the unitary transformation can also be used in other approaches, 
for instance in high-order series expansions by orthogonal transformations ||86l . to capture 
resonant behavior. Not the life time will be accessible directly by a series, but the effective 
Hamiltonian. Then a subsequent analysis by variational methods (as presented in the present 
work) or by diagrammatic approaches has to be used. 

In the second part of this paper we illustrated the theoretical deliberations for the 
asymmetric antiferromagnetic spin- 1/2 Heisenberg ladder. This model shows spontaneous 
triplon decay into two-triplon scattering states. The strength of this decay is controlled by the 
frustrating diagonal interaction J^g- A more comprehensive study of this particular model 
is left to ongoing research. 

We used the generator Fg$(l) which only isolates the ground state and an additional 
Lanczos tridiagonalization in a variational subspace which consisted of states containing up 
to three triplons. We showed that in this way the resonance behavior of the decaying triplon 
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can be described explicitly. The continuous unitary transformation was realized in a self- 
similar way. 

In conclusion, we extended the range of applications of continuous unitary 
transformation to systems which exhibit unstable quasiparticles. We expect that an analogous 
extension can also be implemented for other unitary transformations. 
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Appendix A. Ordering of the generator FpcCO 

The generator FpcQ) sorts the eigenvalues E t in ascending order of the particle number of 
the corresponding eigenvectors such that > qj => E t > Ej. Particularly, this implies that the 
vacuum state |0) of the effective Hamiltonian represents the ground state of the system. In the 
following, we derive this statement. 

In an eigenbasis of the operator Q, which counts the number of (quasi)particles of a given 
state, the generator Fp C (/) is given by 

r i ,j(Q = sgn(q i -q J )<H i , J (r), (A.l) 

where qt and qj are eigenvalues of the operator Q. In general the eigenspace for a given 
number of (quasi)particles qj has a large dimension which is infinite for infinite system size, 
but for the purpose of the present derivation we stick to finite dimensional Hilbert spaces. We 
use the convention that ( Hij(l) stands not only for a single matrix element but for the whole 
submatrix of the Hamiltonian 'H which connects the eigenspace belonging to the eigenvalue 
q t with the eigenspace belonging to the eigenvalue qj. Therefore 'H^ff) is given by a matrix 
with the dimension D, x Dj, where D, is the dimension of the eigenspace Q = q t . 

Using the generator (|A.1I) the general flow equation (12.11) yields the matrix equation 

dfHijiD = -sgn (q t - qj ) (W y (Z)7/y(Z) ~ HjiWjjiD) 

+ £ (sgn (q t - q k ) + sgn ( qj - q k )) <%(/)<H* >7 (0. (A>2) 

Since the effective model will be block diagonal, all off-diagonal matrices < H,j(/) with i + j 
have to vanish for / — » oo. Hence for large / the equation (|A.2I) is dominated by the first term 
on the right hand side where the off-diagonal matrices only appear linearly. So for large / the 
asymptotic behavior of (|A.2I) is given by 

d{H u m = -sgn (qt - q j) (HjHjd) ~ HjiD^j) + O (^(O) . (A.3) 
Note that within this approximation d{H^(T) = Vi, so that we can neglect the /-dependence 
of HiJS) an d *Hjj(l). Without loss of generality we assume in the following that q t > qj. Then 
(TA3T) yields 

dfH u (l) = - ('Hu'Hu(l) - MijWijj) + O . (A.4) 
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The matrix "K, , and the matrix <Hjj are Hermitian, thus unitary transformations 1Y, and 11 j 
exist which diagonalize < H i<i and < H h j, respectively. We will denote this diagonal matrices by 
T>{ := HjHjKi and £> ; := K]<Hjj1lj. By multiplying dA4l) from the left by lij and from 
the right by 1/ ; one obtains 

d,nj(l) = - (Di<Hi,j(D - n,j(D£>j) + O (<%/)) , (A.5) 

where := K] 1(^(1)11 j. According to (|A.5I) the (n, m) matrix element of "H;,;© satisfies 

* * (A.6) 

= -(fo) -fe) 

\\ V J 'm,m) \ J >n,m 

in linear order in the non-diagonal matrices. Since *Hij(l) vanishes for / —* oo, < H i j{l) must 
vanish as well. Therefore, for large / the inequality 

(£>,•) -(©;) >0 (A.7) 

must be fulfilled for all n, m for which the matrix elements CHij(l)) n ,m & are non-zero. 
This implies that all eigenvalues of Di must be larger than the eigenvalues of Dj. Thus, 
the eigenvalues are sorted in ascending order of the particle number of the corresponding 
eigenvectors, as asserted above. 

Note that this ordering does not need to occur, if the corresponding eigenvectors are not 
connected to each other by a finite matrix element of the Hamiltonian for / large, but not 
infinite. If H^fif) = for all / or for I > l < oo the argument to derive (IA.7I) from (IA.6I) does 
not hold. For example, this is the case when the system exhibits symmetries which prevent 
certain subspaces to be linked, as we see for the symmetric spin ladder. 



Appendix B. Transformation of subspaces 

Appendix B .1 . Ground state 

In section [241 we argue that all generators Fp C (l), Fg S (l), F m gsQ), and Fg S ip(/)) considered 
so far transform the vacuum state |0(/)) in the same way if the flow equation is solved exactly. 
Here we prove this statement. 

Previously, we defined the /-dependent Hamiltonian by H(l) := U\l)HU(l). 
Alternatively, we can keep the operators constant but make the states /-dependent. This is in 
complete analogy to passing from the Heisenberg picture to the Schrodinger picture. Hence, 
the /-dependence of the vacuum state is given by |0(/)) = U(l) |0) and the generator is given 
by F(l) = -U\l) (d t U(l)). Thus, for the derivative of |0(/)> it follows 

d l \0(l)) = d I U(D |0> 

= U{T)U\l) (<9,[/(Q)|0> 

v v ' (is.!) 

=-F(l) 

= -U(l)F(l)\0). 
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Introducing a basis {|z)} yields 

d, |0(0) = - J] U(l) \i) <£|_F(0|0> . 

The key observation is that for all considered generators the matrix element F ifi (l) is the same, 
namely 

(H i0 (l) for*>0 

F '» m = {o' fori = o- <B ' 3) 

Applying (lB3l) to (lR2l) yields 

3/|0(0> = -£t/(OIO<*l#(0 10) 

( \ (B.4) 

= - YjU(D\i)(i\H(D\0) +t/(0|0)<0|^(0|0). 

Shifting the /-dependency to the vacuum state and using the equality 17(1)17^(1) = 1 provides 
us with 

dt |0(/)) = -H |0(0) + 10(0) (0(01 H |0(0) 

r n (B.5) 

= [p (D,h] |0(0) 

with the /-dependent projector Pq(0 = |0(/)) (0(/)|. According to (IB.5I) the derivative of 
|0(/)) only depends on |0(/)) itself and the initial Hamiltonian H. Therefore, the considered 
generators all transform the vacuum state |0(/)) in the same way. The essential point of the 
proof is that for all considered generators the matrix elements F !i o(0 are defined identically 
by (IB.3I) . Note, however, that the statement, that all generators treat |0) alike, does no longer 
hold if approximations (truncations) are introduced. 

Appendix B.2. One-particle space 

The proof presented in the previous subsection can be generalized. Since the action of the 
generator i 7 pc(/) and the generator Fg S ^ p(0 is also the same on the one-particle subspace, one 
can prove that they also transform all one-particle states in the same way. In the following we 
characterize the states by their number of (quasi)particles, so it is useful to use an eigenbasis 
{|/)} of the (quasi)particle number operator Q. The number of (quasi)particles of a state \i) is 
denoted by g,. Consider the derivative of all states with at most one particle 

= -J]u(i)F(im 

i:?f<l 

= -YjU(d \j)<j\Fm}. 

For both generators the matrix elements Fjj(l) with q, < 1 are given by 

Fjj(l) = sgn(qj-q,)Hjj(l) (B.7) 
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according to (12.31) and (I2.20I ). Hence we have 

dij] \m = -Y J U(D\j)(j\H(D\i) 

i:qi<\ m<\,j:qj>l ^ ^ 

-Y^sgn(qj- qi )U(r)\j)(j\H(D\i) 

i:qi<\,j:qj<\ 

To the first part on the right hand side of (|B.8I) we add all missing contributions with q } < 1. 
Hence we arrive at 

i:qi<\ i: qi <\ i:qi<\,j;qj<\ ^ 

-Y J Sgn(q J -q l )U(D\j)(j\H(D\i). 

i:qi<l,j:qj<l 

Just as in the previous subsection we shift the /-dependence from the Hamiltonian H(l) to the 
states 

diYj l*'(0> = ~Tj H l*'(0> + Z ^ ,(/)) w H |/(/)) 

(':oi<l i:<?i<l i:a;<l,/:<7;<l 

* (B.10) 
-2,sgn(^-^)|j(/))0'(/)l^l^)). 

i:qi<l,j:qj<l 

It follows that the transformation of the subspace {|z)} with q t < 1 is independent from all 
other states {|z')} with g, > 1. The transformation only depends on the initial Hamiltonian 
H. Therefore, the generator i*pc(0 and the generator Fg S ip(/) transform the one-particle 
subspace in the same way. Note that this proof is not restricted to the case < 1 and can 
easily be adapted to the case qi < n e N. The choice (12.201) or (IB.7I) has to be adapted 
accordingly, i.e., we have pass from Fg S ip(Z) to Fg SjW p(/) with 

n N 
p=0 i>/? 



Appendix C. Lanczos tridiagonalization 

To determine the spectral density (I4.7I) one has to calculate the coefficients ai(K) and 
of the continued fraction representation of the Green function (I4.8I) . This can be done by 
the Lanczos recursion scheme, for which we restrict our calculations to the subspace (I4.9qI) - 
(|4.9cl) . Within this subspace the recursion (Lanczos tridiagonalization) 



l«Ao> = 


\K,z) 






(C.lfl) 


l<Ai> = 


/r/res 
reff 


- a (K)) 




(C.lfc) 


Wi) = 


/r/res 
l^eff 




-^iW 2 W 


(C.lc) 


l<A 3 ) = 


/r/res 
l^eff 


-a 2 (Z))|<A 2 > 


- W) 2 l«Ai) 


(C.1J) 
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Figure CI. Continued fraction coefficients cij(K) and bj(K) for x = 0.5, y = 0.1 and 
momentum K = 7r/10. The restriction of the considered subspace becomes conspicuous at 
i ~ 115. 



with 



am = 

bi(Kf = 
bo(K) 



<(A,_i|iA/_i) 




for? = 0,1,2,... 



for i = 1, 2, 3, 



(C.le) 

(CI/) 
(Clg) 



generates a set of orthogonal states The action of the restricted Hamiltonian (14.101) in the 
three-particles subspace (|4.9<ab - (|4.9cb is given in section Appendix D In the generated basis 
{\iffi}} the matrix of the restricted effective Hamiltonian is tridiagonal, where the a,(X) 
are the diagonal matrix elements and the bi(K) are the elements on the second diagonal. All 
other matrix elements are zero. 

Figure ICTI shows the results for the coefficients ai(K) and bj(K) for x = 0.5, y = 0.1 
and K = n/10. First it appears that both coefficients a,(/Q and b t (K) (i ^ 115) converge to 
fixed values a^K) and boo(K) as it should be for a bounded and gapless spectral density of 
an infinitely large system [1841 . But for i ^ 115 both coefficients start to change their values 
again noticeably. This is a consequence of the fact that we had to restrict the relative distances 
d and d\ + d 2 to a maximum of 119 rungs (cf. (|4.9al )- (|4.9cl )) in our numerical calculations. 
Therefore, we only use the first 100 coefficients and terminate the continued fraction at i = 100 
as described in the following subsection. 



Appendix C.l. Termination 

The spectral density S(K,to) at fixed K as obtained by a finite continued fraction of the 
Green function (|4.8I) has poles at the zeros of the denominator. Thus, the spectral density 
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S(K,to) is a collection of 5-peaks. One standard approach to obtain a continuous density is 
to introduce a slight broadening of S (K, to) via to — > to + i<5 with a small real number 6. This 
procedure corresponds to smearing out 5-peaks as Lorentzian functions of width 6. The caveat 
is that also all truly sharp features such as band edges or van Hove sigularities are smeared 
out. However, a notably improved resolution of S (K, to) can be achieved by introducing an 
appropriate termination of the continued fraction. 

If we want to evaluate the continued fraction for the infinite large system, we have to 
stop our recursion before the finiteness of the considered subspace (I4.9al) - (l4.9cl) becomes 
conspicuous. Therefore, we compute the average value of aj(K) and bj(K) for i = 80 ... 100 
(cf. figure ICTT) to obtain a good approximation for the limits aoo{K) and bco(K). From these 
limits the upper and lower boundaries E U (K) and E^(K) of the spectral density S (K, to) are 
determined via [ 841 



E U (K) = aoo (K) + Ib^K) (C.2a) 

E X (K) = aUK) - 2Z?co<X) • (C.2b) 

The existence of an upper boundary of the spectral density is a consequence of our restriction 
to a subspace which contains three quasiparticles at maximum. Finally, we use the directly 
calculated coefficients ai(K) and b i+ i(K) for i = 0...99 and subsequently the square root 
terminator defined by the approximate limits a^K) and b^K). By doing this, we assume 
that all following coefficients at(K) and b i+ i(K) with i > 100 are constant. The square root 
terminator T(K) is given by 

T = — r (to- a^- V^D) for to > E u (C.3a) 
v ; 

T = -^r- (co-a^-i Vd) for & < to < E u (C3b) 
Zbl v I 

T = ^—(co-a oa + V^D) for to < E { (C.3c) 



with 



D = 4bi - (to - aS , (C.3J) 



where we suppressed all K dependences for the sake of simplicity. The last considered 
coefficient b l00 (K) in (14.81 ) is multiplied by the appropriate terminator T(K). The imaginary 
part of the resulting expression yields the continuous part of the spectral density S(K, to). In 
this way we can reliably approximate the thermodynamic limit of the spectral density S (K, to) 
by calculations in a finite subspace. 



Appendix D. Analysis of the effective model 

Here we present details of the analysis of the effective Hamiltonians for the asymmetric ladder 
generated by CUTs. More details are given in |[87l HH [881 for the special case of a particle 
conserving effective Hamiltonian. 



Adapted CUT to treat systems with quasiparticles of finite lifetime 33 



HQ= oo): 




(a) (b) (c) 



Figure Dl. (Color online) Considered subspaces to analyze the effective Hamiltonian obtained 
by the generator Fgs(l). The generator Fgs(l) only isolates the part H^(l). Colored blocks 
illustrate interactions which are included. Grey blocks illustrate neglected interactions. Panel 
(a) describes an analysis within the one-triplon subspace. Panel (b) describes an analysis 
within the one- and two-triplon subspace. Panel (c) describes an analysis within the one-, two- 
and three-triplon subspace, see ( |4.9al >- d4~9cb 

The generators Fpc(/) and Fg S jp(/) isolate the one-particle subspace from all other 
subspaces (cf. figure [2k and figure [2f). Therefore, the one-particle eigenvalues can be 
calculated without considering states with a higher particle number. 

The effective Hamiltonians obtained by the generator Fg S (l) still contain interactions 
between the one-particle subspace and other subspaces (cf. figure [2b). Consequently, a 
diagonalization in the one-particle subspace only gives an approximation for the eigenvalues 
of the effective Hamiltonian, namely an upper bound for the eigenvalues if the ground state 
energy is sufficiently well described. The results for the eigenvalues can be improved by 
considering higher particle subspaces as well (see figure ID1I) . In this paper we consider 
subspaces which consist of states which contain up to three triplons. 

The Fourier transformed one-, two- and three-particle states are given by (|4.9al )- (|4.9cl) . 
In any practical calculation the relative distances must be truncated to make the subspace 
finite. Note that due to the hard-core algebra it is not possible that two particles occupy the 
same rung. Below the action of the various parts of the effective Hamiltonians are given. 

Appendix D.l. H\ 

The action of the operator H\ on the one-triplon state \K, a) is given by 

H\\K,a)= Yj^ Kr c1C\K,a') (D.l) 

r,a' 

with 

c a u f r := (R,a'\H{\R + r,a) (D.2) 
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and R e Z. Note that due to the translational symmetry only the relative distance r between 
the states \R, a') and \R + r, a) occurs in the coefficient c"£. In the special case of a SU(2) 
symmetric Hamiltonian the coefficients c n '" obey the relation c n '. a = S a ',a c n-,r- The used 
truncation scheme (cf. section 13711) causes c n '" = for |r| > J 2 - All other coefficients c which 
appear in the following are affected by the truncation scheme in an analogous way. 
The action of the operator H\ on the two-triplon state \K, a) \d,/3) is given by 

H\ \K, a) \d,j3) = ^ cfi^e** \K, a') \r + d,/3) (D.3a) 

r>-d 
a' 

a' 

+ E ^ a) |_ (r - J) (D3c) 



r>d 



The action of the operator H\ on the three-triplon state \K, a) \d\,/3) \d%, y) is given by 
H\ \K,a) \d u /3) \d 2 ,y) = E cf^e** \K,a') \r + d u /3) \d 2 ,y) (DAa) 



r>-di 
a' 



h+d2)<r<-d\ 
a' 

+ E c n" eiK * \ K >P> r> I- (' + * + 40 , <*') (D.4c) 

-<-(t/i +d 2 ) 
a' 

—di<r<d\ 
P 

+ E ^i^ 5 ^ a) 1 Ji + r) |_ (r + J2) (DAf) 

r<-di 
P 

+ E C S eiAr§ l^«>l^i^>|-(r-rf 2 ),y> (D.4g) 

r' 

+ E C nl QlKh {K > a) l_ (r " dl ~ dl) ' y,) |r " (D - 4/z) 

\l<r<d\ +d2 
7' 

+ E c ^ e ^ I*. /> l r - * - rf 2, |rfi,j8> (D.4i) 



r>di+do 

Y 
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Appendix D. 2. H\ 

The action of the operator H\ on the one-triplon state \K, a) is given by 

H\ \K, a)=J] 4%<rtrQ \K, a') \d',/3') (D.5) 



r,d' 

with 

v on zj2 



The action of the operator H\ on the two-triplon state \K, a) \d,fi) is given by 

Hj \k, a) \d,p) = J] 4it;> lK[ ^ ) i^. <*') k*'. /> \ r+d ~ d '^) (°- 7 «) 



r>-d+d' 
d',a',y' 



d<r<-d+<i' 
d',a',y' 

+ 2 4';^^) |*,0> I- (r + d), a') \d', y') (D.7c) 

r<-d 
d',a',y' 

+ Z ^iS^*^ I*. «> l _r + <W \d'> J) (D.ld) 

r<d 
d'fi',y' 

+ E 4^^^) | *,/0 |r - J, a) |-r + J + J', /> (D.7«) 

+ 2 cfi^^^) l^,/i'> y'> \r-(d + d') , a) (DJf) 



r>d+d' 



Appendix D. 3. H\ 

The action of the operator H\ on the two-triplon state \K, a) \d,/3) is given by 

H\ \K, a) \d,/3) = 2 c^e K ^ \K, a') (D.8) 



r,a' 

with 

y,afi ._ /p /, rrl 



c Tw := a '' H 2 I* + «> I* + r + d,/3> • (D.9) 
The action of the operator H\ on the three-triplon state \K, a) \d\,f}) \d 2 , y) is given by 

H l 2 \K,a)\d u p)\d 2 ,y) = £ c^e a'> |r + ^ + J 2 ,y> (D.lOa) 



r>-(rf,+J 2 ) 
a' 



c iCi e \ K > I" ( r + d i+ d 2) , a') (D.lOfc) 



r<-(rfi+rf 2 ) 
a' 
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, . (3r+d l +2d 2 \ 



r<d\ 



. I 3r+rf[ +2d 2 \ 

+ J]^/ K \^)\K,/3')\r-d u a) 

r>d\ 
P' 

( 3r+d l+ 2d 2 \ 

r>-d\ 
a' 

, ( 3 r +d l+ 2d 2 \ 
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(D.lOc) 
(D.lOd) 
(D.lOe) 
(D.IO/) 



r<-d\ 

a' 



Appendix DA. H\ 

The action of the operator H\ on the one-triplon state \K, a) is given by 



2d' +d' 



Hl\K,a)=J]c^;^JT 3 )\K,a')\d[,p)\d' 2 ,Y) 



r,d[,d' 2 



with 



4£££ : = <*> a '\ < R + d 'i>P\ < R + < + 4. r'l #i I* + r, a) 



(D.ll) 



(D.12) 



Appendix D. 5. H\ 

The action of the operator on the three-triplon state \K, a) \di,/3) \d 2 , y) is given by 

H\ \K, a) \d u /3) \d 2 , r > = 2 Ca^' a '> (D - 13) 

with 

c f 3 j£fe := «'l H l \R + r,a)\R + r + d u fi \R + r + d x + d 2 , y) . (D.14) 

Appendix D. 6. 

The action of the operator H\ on the two-triplon state \K, a) \d,/3) is given by 

H\ \K, a) \d,fi) = J] 4SrP^^ \K, a') \d',/3') (D.15) 



r,d' 
a'£' 



with 



4&$4 '■= o'l + d '>P'\ H l \R + r,a)\R + r + d,{3) . (D. 16) 

The action of the operator H\ on the three-triplon state \K, a) \d\,fi) \d 2 , y) is given by 



Hl\K,a)\d u /3)\d 2 ,y) = J]c^e 

r>-(d\ +d 2 )+d[ 
d[,a'J3' 



2r+d { -d' 



1) 

l\K,a')\d[,fi')\r + d l +d 2 -d[,y) (D.ll a) 
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I 2r+d{ -d'j 



+ V c a ' ,fl '' a 

/ i 22\r,d\ ,d\ 
-{d\ +d 2 )<r<-(d\ +d 2 )+d' 1 
d[,a'fi' 
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fi Q \ 3 l\K,a')\r + d l +d 2 ,y)\-(r + d l +d 2 ) + d[,p')(D.nb) 



r<-(^i+d2) 
d' v a'J3' 

T /_j C 22;r,rf 1 +^,d 1 +d 2 
-iil<r<rf2 



2r+d 2 -<£ 



i|#,aO|r + Ji,/3>|-r + 4,y'> 



+ 



/ i ^22;r,, 



d\ +d'd[+d 2 



r<-d[ 
d' 2 ,a',y' 



+ 



\K,P) |- (r + d{) , a') |Ji + d' 2 , y) 

I 2r+d 2 -d' 2 \ 

'A 

l\K,a)\-r + d u p)\d 2 ,y') 



d' 2 ,a',y' 
+ 



22;r,d' 2 , 



r<d\ 

d' 2 ,P',r 



2 4 2 



+ 



22;r,d' 2 ,. 



2r+do -dt 



di<r<di+d' 2 

d' 2 fi',y 



V 3 /|tf, ) S'>|r-<* 1 ,a>|-r + <* 1 +4, r '> 



+ 



22;r,d' 2 , 



.l 2r + d 2 -d' 2 \ 

■ge \ 3 i|^ J S')|4,r'>|r-(J 1 +4),or) 



r>di+d' ? 



(D.17c) 



(D.17J) 



(D.17e) 



(D.17/) 



(D.17g) 



(D.17/i) 



(D.17/) 



Appendix D. 7. H\ 

The action of the operator H\ on the two-triplon state \K, a) \d,/3) is given by 

/ Id-Ad' -2d' \ 

Hl\K,a)\d,p) = J]c a ^JJ^ ' 6 i|^a')K^')l4, r ') 



with 

ot'$',y\<xJS 
t 32;r,d' v d' 2 ,d 



(R, a'\ (R + d\,p\ (R + d\ + d' 2 , y'\ H\ \R + r,a)\R + r + d,p) 



Appendix D. 8. Hj 



(D.18) 



(D.19) 



The action of the operator Hj on the three-triplon state \K, a) \d\,j3) \d 2 , y) is given by 

^— , , , tA *d l +2d 2 -3d' \ 

H\ \K, a) \d u P> \d 2 , y) = Yj c &j5Zh e V+ 6 ' \ K - a '> l d '^'> ( D - 2 °) 



r,d' 
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Table Dl. States of the 5 = 1 and m — subspace. 

|5 = l,m = 0> 1 |z> 

|S = l,m = 0> 2 ^(|*,y>-|y,*>) 

\S = l,m = 0) 3a -^(\z,x,x) + \z,y,y) + \x,z,x) + \y,z,y)J 

~ yji ( 2 l z ' z ' z > ~ I x ' x ' z > ~ tV' y. z> ) 
|S = l,m = 0) 3i ±(\z,x,x) + \z,y,y)-\x,z,x)-\y,z,y)j 
|S = l,m = 0> 3e ^(|jc,je,z> + |y,y,z> + |z,z,z>) 



with 

4Saa := a '' (i? + #3 I* + r, a> |* + r + J 1? /?> |* + r + d x + d 2 , y) . (D.21) 



Appendix D. 9. 

The action of the operator H\ on the three-triplon state \K, a) \d\,f3) \d 2 , y) is given by 

i X + 2 ( rf i- / iM d 2-4)l 

Hi \K, a) \d u P)\d2,j) = J] 4^AA e 1 i|^a')K,y8')l4,r') (D.22) 

a',/?',/ 

with 

: = («. «'l + 4>fi\ ( R + < + 4. r'l flf 

^ " (D.23) 
x |7? + r, a) \R + r + d u /3) \R + r + d\ + d 2 , y) . 



Appendix D. 10. S = 1, m = subspace 

Due to the SU(2) it is possible to reduce the computational effort. The one-, two- and three- 
triplon states with 5 = 1 and m = are listed in table ID 1 1 Since they are independent from 
the total momentum K and the relative distances d, d\ and d 2 we omit the dependence on these 
parameters. 
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